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Abstract 



We determine the thermodynamic stability conditions for near-extreme 
rotating D3, M5, and M2-branes with multiple angular momenta. Critical 
exponents near the boundary of stability are discussed and compared with 
a naive field theory model. From a partially numerical computation we 
conclude that outside the boundary of stability, the angular momentum 
density tends to become spatially inhomogeneous. 

Periodic Euclidean spinning brane solutions have been studied as mod- 
els of QCD. We explain how supersymmetry is restored in the world- 
volume field theory in a limit where spin becomes large compared to 
total energy. We discuss the hierarchy of energy scales that develops as 
this limit is approached. 
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I. INTRODUCTION 



Mathematically, thermodynamic stability is equated with the subadditivity of the 
entropy function. More precisely, suppose the entropy S is known in the microcanonical 
ensemble as a function of the other extensive thermodynamic variables (like energy and 
charge), which we will collectively denote Xt. Subadditivity of the function S({xi}) is 
the condition 

S({\x t + (1 - A)5,}) > XS({xi}) + (1 - X)S({x i }) . (1) 

If the inequality went the other way, then the system could gain entropy by dividing 
into two parts, one with a fraction A of the energy, charge, etc., and the other with a 
fraction 1 — A. When the system satisfies (1), no such process is allowed by the Second 
Law of Thermodynamics, and it is in this sense that the system is thermodynamically 
stable. 

For black holes, perhaps a better way to think of (1) (though a more heuristic one) 
is as follows: (1) is satisfied for all Xi and X{ in a sufficiently small neighborhood of a 
given point in phase space if and only if it is possible for a black hole at that point 
to exist in thermal equilibrium with a heat bath much larger than itself. A simple 
example would be the Reissner-Nordstrom black hole in four- dimensional Einstein- 
Maxwell theory with no electrically charged particles. Near extremality, the specific 
heat is positive, but far from extremality the specific heat is negative. If one starts 
with the black hole in a state with positive specific heat and a temperature close to the 
temperature of some thermal gas of gravitons and photons, then the black hole will 
equilibrate with the thermal gas. On the other hand, if the black hole was in a state 
of negative specific heat, then even if its temperature is arbitrarily close to that of the 
thermal gas, the black hole will either emit faster than it absorbs or absorb faster than 
it emits, and correspondingly get smaller and hotter or bigger and colder. There is no 
bound on the runaway phenomenon except for the black hole to make the transition to 
positive specific heat or else to absorb most of the matter in the universe (assuming the 
universe is finite). That is what one means by instability of the initial configuration. 1 

For spinning D- or M-branes, the heat bath line of reasoning is questionable be- 
cause it is not clear whether one can construct an appropriate heat bath. One could 
perform the following gedanken experiment. Start with an array of nearly identical 
parallel spinning branes distributed through space, and ask whether, as they trade 



In this discussion we have assumed that the Jean's instability of the thermal gas is 
insignificant — that is, it enters into the problem with a much longer time-scale than the 
time scale of equilibration between the black hole and the thermal gas. 
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Hawking radiation among themselves, there is a process of equilibration or disequili- 
bration among them. This is the standard trick of constructing an ensemble out of 
replicas of a given system. The branes should equilibrate precisely when (1) is sat- 
isfied. However, there is a difficulty: one has to estimate whether the time scale of 
equilibration or disequilibration is faster or slower than the time scale on which the 
gravitational forces between the replicas appreciably distort the array. Near the phase 
transition of [1], one expects the equilibration time scale to diverge, so the gedanken 
experiment will become impractical (even by the standards of gedanken experiments) 
close enough to the transition. 

As a first step we can think of a single cluster of near-extreme coincident branes 
in a decoupling limit of the theory, where Hawking radiation to the asymptotically 
flat region of spacetime is very slow. Once we have analyzed the criterion (1) and 
established the existence of unstable solutions (even in this near-extreme decoupling 
limit), we can ask what form the instability takes. Because the branes have spatial 
extent, perhaps the most natural assumption is that the solutions which spin too 
fast are unstable against distributing their angular momentum unevenly across their 
worldvolume, somewhat analogous to the Gregory-Laflamme transition [2]. We will 
have more to say on this possibility in section III. 

The work of [1] illustrates that the thermodynamics of the D3-brane world-volume 
theory, as read off from the Bekenstein-Hawking formula from the supergravity geom- 
etry, is stable only on a restricted region of phase space. That is, (1) is guaranteed 
to be satisfied only if the Xi and the Xi lie within the specified region. Our first goal, 
which we pursue in section II A, is to extend the stability analysis of [1] from one inde- 
pendent angular momentum to three, and to the case of near-extreme M5-branes and 
M2-branes with two and four angular momenta, respectively. (The maximum number 
of independent angular momenta is the rank of the rotation group perpendicular to 
the brane). 

We will restrict ourselves to the near-extreme limit of spinning branes. For D3- 
branes, the string theory account of the microstates of the brane in terms of strongly 
coupled Af = 4 super- Yang-Mills theory is thought to apply in this limit. The near- 
extremal spinning geometries have a near-horizon limit which is asymptotic to AdS 5 x 
S 5 , so if we employ the philosophy of the AdS/CFT correspondence [3-5], we would 
say that the supergravity thermodynamics translates directly into a characterization of 
thermal states of the gauge theory. Instability is more of a surprise in this context, since 
statistical mechanical derivations of thermodynamics as a rule imply subadditivity of 
the entropy. The difficulty of treating the strongly coupled Af = 4 gauge theory at 
finite temperature prevents us from deriving or disproving the instability from first 
principles. In section II C we do however extend the simplified field theory treatment 
introduced in [1] to the case of several angular momenta, and we show that the results 
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are qualitatively similar to the supergravity predictions. 

In section II B we discuss the supergravity predictions for critical exponents at the 
boundary of stability The maximal number of fractional critical exponents at a given 
point on the boundary is the number of non-zero angular momenta. We show in 
section II C that similar critical exponents arise from the naive regulated field theory 
treatment for the D3-brane. Our analysis is not detailed on this point. 

In section IV we analyze the supergravity solutions Wick-rotated to Euclidean time. 
We observe that the angular momentum amounts in field theory to a twist imposed on 
the -R-charged fields in the field theory: the field theory partition function is altered 
from tre _/3H to tr Re~^ H where R is an element of the i?-symmetry group. It describes 
the angle through which the brane spins during one period of the Euclidean time. In the 
limit discussed in [6] where angular momentum becomes large with energy held fixed, 
R — > (— 1) F , so we wind up with almost supersymmetric boundary conditions. For M5- 
branes compactified down to four dimensions, this leads to Af = 4 super- Yang-Mills 
theory with supersymmetry slightly broken by boundary conditions in the Kaluza- 
Klein reduction. We review how the near-supersymmetry affects the proposals to use 
spinning branes as holographic approaches to QCD, and we note a puzzle regarding 
the energy scale of supersymmetry breaking and confinement. 

II. THERMODYNAMICS OF SPINNING BRANES 
A. Supergravity analysis 

The supergravity solutions for spinning branes were constructed in [7]. 2 The metric 
for the general spinning D3-brane and M5-brane can also be found in [9,10] and [11], 
respectively. We are interested primarily in the thermodynamics, and this follows 
directly from the general formulas of [8] . 

Some parameters that enter into the thermodynamics of D3-, M2-, and M5-branes 
are shown in the following table: 



2 Unfortunately, the spinning M-brane solutions quoted in [7] have a number of typographical 
errors. However, these general spinning M- as well as D-branes can be obtained in a straight- 
forward manner from [8]. There the general D-dimensional rotating two-charge black hole 
solutions (of toroidally compactified string theory) were given. The corresponding spinning 
(D3, M5, M2)-branes are obtained by turning off one of the two charges of the (D = 7, 6, 9- 
dimensional) black hole, and then lifting the solution in a straightforward way to ten and 
eleven dimensions, respectively. 
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brane p D T p a p c p 

D3: 3 7 ^ 2 -L 

/ 7T \V3 1 I 2 ) 

M5: 5 6 — 3 



M2: 2 9 



V2kV 96tt 3 
'2tt 2 V /3 3 1 



/ 2 3^ 



In (2), k is the ten-dimensional gravitational coupling for the D3-brane, and the eleven- 
dimensional gravitational coupling for the M2-brane and the M5-brane. Using formulas 
(14) and (15) of [8], one finds the following expressions for the mass, charge, and spins 
of the branes: 



M 



16.£ (D - 3) ( C ° Sh2a+ ^3 



(3) 



Q = °~ 20 — (D - 3) cosh a sinh a = NT P 
167rG£)_2 

J i = Tc n 2 Vi cosh a 

s = ^D-2To ^ cogh a _ 

167tGd_ 2 

Our variables are defined in terms of those of [8] by 

rt 3 = 2m m = ^ y H = -. (4) 

The location of the horizon, y# , is the largest positive root of a cubic or quartic equation 
(obtained from (16) of [8]): 

D3: (y 2 + yf) (y 2 + y 2 2 ) (y 2 + y 2 3 ) - y 2 = 

M5: (y 2 + yj)(y 2 + y 2 2 )-y = (5) 

M2: (y 2 + yl){y 2 + y 2 2 )(y 2 + y 2 )(y 2 + y\) - y 2 = 

It is useful to introduce one more combination of the parameters: 

R D ~ 3 = r^ 3 cosh a sinh a . (6) 

The constants a p and c p were defined so that 
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D-2 



N a p2K 2 



w forD3 

^ forM2andM5. 



(7) 



We can now define E = M — Q and rewrite 



E 



r 



D-3 



(D — 3) ( cosh 2 a — cosh a sinh a H 

N a rV Rp +1 \RJ y J \ D-3 



c p D-l fr_o\ D ~ 3 
RJ 

D-2 



R P+l 2 



Ji = 



s = 



Ji c p fr 
-— : — — = -^-2vi cosh a — 
N a ?V Rp y \R 

——: — — = -^-4:TTy H cosh a ( ) 
N a ?V Rp y \RJ 



D-2 



(8) 



(D-l)/2 



The approximate expressions are for the near-extreme limit, a — > oo with i? held 
fixed. Scaling out the power iV ap \/ from E 1 , J i; and 5 allows us to write all subsequent 
expression in terms of e, ji, and s alone. That all powers of V can be absorbed in this 
way is dictated by the absence of any scale in the theory; that all powers of N can also 
be absorbed seems to be telling us that the effective number of degrees of freedom is 
on the order N ap , as first conjectured in [12]. We can form dimensionless ratios 



J 



eP 
p+i 



eP 

J± 

s 



= c 



2 2 p +1 (2tt)p + 1 
Cp (D-1)p 

22p+l 



p (D-1)p 

}_Vj_ 
2ny H ' 



p+i 



-Vh 



p+i 
Vi 



(9) 



1. Grand- canonical ensemble 

Now we wish to ask in what region of phase space s is a sub-additive function of the 
energy density e and the angular momentum densities ji. We regard the ji as thermo- 
dynamic variables rather than fixed parameters, and the requirement of sub-additivity 
means that Legendre transforms with respect to the ji as well as with respect to e are 
invertible. Thus, in the sense described in section III A of [13], we are investigating 
stability in the grand canonical ensemble. Stability in the canonical ensemble, where 
only the Legendre transform with respect to e is required to be invertible, amounts to 
positivity of the specific heat, with the ji regarded as fixed parameters. We will treat 



6 



stability in the canonical ensemble in section II A 2. We will exhibit evidence in sec- 
tion III that the grand-canonical ensemble is the relevant one when the world-volume 
is large. All of our actual calculations will be done in the microcanonical ensemble, 
that is writing s = s(e,ji). 

The function yu = yuiVi) is complicated to express in closed form, since it is the 
solution of a cubic or a quartic equation, cf. (5). Using (9), we can write 

. = (^j) , (10) 

where 7 P+1 = c p 2 2p+1 /(D — and the arguments yi of yn have been written out 
explicitly in terms of the ji. Fortunately, there is a more convenient way to judge 
the subadditivity of (10) than brute force. For instance, one can require that e is 
concave up as a function of s and the jV Or, most usefully, one can show that the 
domain of stability is the maximal region including the points where all ji = where 
the temperature T = de/ds and the angular velocities f2j = de/dji are one-to-one 
functions of s and the This is the criterion we will actually use. 

By differentiating the explicit expression (10) with respect to ji, one can show that 

It is useful to define the dimensionless ratios 

Xi = J*. = 2tt^ . (12) 
Vh s 

Because of the scaling which follows from conformal invariance, any dimensionless 
thermodynamic quantity (like e p /s p+1 ) can be expressed in terms of the Xi. This is 
true in particular of Qi/T: we have 

^L = -^W where W = 2 - . (13) 

oyi 1+a^i 3 — D + 2^-1 Xj 



The sum over j in W is over all the independent angular momenta, of which there are 
\_(D — 1)/2J. The most straightforward way to prove (13) is to use the relations 

^ = 11(1 + *?) 

i 

(3- )f = E^ (14) 

ob^ _ dm. _ dyH_ 

Xi yi y H 
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The first of these is a rephrasing of (5); the second follows from the first by logarith- 
mic differentiation; and the third follows from (12) also by logarithmic differentiation. 
Again because of the various scaling laws, the relations T = de/ds and Jlj = de/dji 
are one-to-one if and only if (13) can be inverted to express the x^ in terms of the 
dyn/dyi. At the boundary of a region of x^-space on which (13) can be inverted, at 



least one of the eigenvalues of the matrix 
can show that W is finite on the entire region o 



d dy H 

dxi dyj 



d dy H 

dxi dyj 



goes to zero or infinity. Provided we 
stability, its boundary will be given 
closest to the origin, X; = 0. Extracting an overall 



by the solution of det 

factor of Vt/L( D_1 )/ 2 J ~[\ k -i, the determinant equation becomes 



D(xi) = det 



d log W 1 - x\ 

Xi ~ "HZ '. o 



dxj 



xt 



= 



(15) 



This determinant has a simple closed form: 

(L>-3)P(-1) + 2P'(-1) 



D(x 

where -P(A) is defined by 



(D — 3)P(1) — 2P'(1) 



LVJ 



(16) 



p(A) = na + a*2) = i + e a £ e n 



;i7) 



=1 ii>t2>--->ii j=l 



Note that D(x{) = 1 when all the Xi vanish. Thus the region of stability is the domain 
including Xi = where D(xi) > 0. If all but one angular momentum is zero, then 



D(x) 



(D-3)-(D- 5)x 2 
(D - 3) + (D - 5)x 2 



(18) 



Which indicates that the region of stability for a single angular momentum is x < 



D-3 



Apparently by coincidence, a v 



D-3 
D-5 



as well. The bound x 2 < a p is in agreement with 



equation (17) of [1] in the case p = 3, and with the bounds derived in [14] for p = 2 
and p = 5. 

It is straightforward to verify that the denominators of W in (13) and D{xj) in (15) 
differ only by a factor — P(l), and that both these denominators are nonzero and finite 
throughout the domain including Xi = where the numerator of D(xj) is positive. 
Thus an equivalent condition of stability to D{xi) > is 

(D - 3)P(-1) + 2P'(-1) = (D - 3) + ± (-l) e (D - 3 - 2£) £ ft x\ > 

1=1 i\>i2-->ii 3=1 

(19) 
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For the three types of branes we are considering, (19) specializes to 



D3: 



2 2 2 2 I 2 2 2 \_ r\ 

X ~* X-^X^X^ _ U 



M5: 



> o 



M2: 



3 - 2(x? + Xa + xl + xl) 



(20) 



■ / 2 2 I 2 2 I 2 2 I 2 2 I 2 2, 2 2\ 
\ yX-^X<2 \~ X-^X^ \~ X -^X ^ — r~ X^X^ ~\~ X^X^ ~\~ X^X^j 



2 2 2 2 r\ 
X-^X^X^X^ ^_ U . 



Our methods can also be used to judge the stability of rotating black strings in five 
and six dimensions arising from intersections of spinning M-branes: M5 _L M5 _L M5 
for five dimensions, and M5 _L M2 for six dimensions. Because the world-volume is 
effectively two-dimensional, one does not expect any phase transition at finite temper- 
ature, on account of the Mermin- Wagner theorem [15]. Indeed, it is straightforward 
to verify by direct computation that the entropy function is subadditive. It is also 
possible to trace through an analysis like the one in section II A, introducing a combi- 
nation of parameters Rf ~ 3 = r ~ 3 cosh on sinh a.% for each constituent family of branes, 
and verifying that (16) is always satisfied. (We note however that (19) is not always 
satisfied for the six-dimensional case: a zero in the denominator of (16) cancels a zero 
in the numerator.) 



For the sake of completeness we also derive the subadditivity constraints for the 
canonical ensemble. This is the case where the angular momenta ( "charges" ) of rotating 
branes are kept fixed, so that subadditivity of the entropy is exactly equivalent to 
positivity of the specific heat at constant j«. As we will demonstrate, it may be difficult 
to keep the j% uniform over the entire world-volume outside of the grand-canonical 
region of stability. However, uniform j% will be our working assumption in this section. 
It is a general result, demonstrated in section III A of [13] and verifiable directly from 
comparing (22) to (20), that the region of stability in the canonical ensemble is larger 
than in the grand canonical ensemble. The canonical stability constraints for spinning 
branes with one angular momentum were first studied by Cai and Soh in [14]. 



Using (8), (9), (10), (13), and (14), it is fairly straightforward to reduce the inequality 

d 2 s/de 2 < to 



2. Canonical Ensemble 



3 3 3 



3 



D3:2 + 55>2-5>f + 10]r 



x 2 x] + l<ox\x\x\ + 2 53 




i=l i=l i<j 
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3 

+3 E x\x A x\ + 5x 1 x 2 x 3 > , 
{i<j}^k 

M5 : 9 + 36(xj + x\) - %x\ + x A 2 ) + 76x 2 x 2 + %{x\x\ + x A x\) + 2lx\x\ > , 

M2 : 9 + 18 E A + 34 E + 60 E z?a^ + 2 E - 3 E a£r* + 6 E afrM 

i=l i<j i<j<k i^j i<j {i<j}^k 

4 4 
_I_QC 2 2 2 2 9 \ 4 4 2 , in \ ~> 2 2 2 4 

-\-i)OX 1 X 2 X 3 X 4: Z x i x j x k' iZ 2^ x i x j x k x l 

{i<j}^k {i<j<k}^e 

+2 E ^x}^ + 6 E ^K^ 2 + 9^444 > , (21) 
{i<j}¥={k<e} {i<j<k}^£ 

(22) 

where the defined in (12). 

A special case a single angular momentum turned on, say, X\ = x ^ 0, these stability 
constraints reduce to: 

D3 : 2 + 5x 2 - x A > , 



M5 
M2 



1 + 4x 2 - x A > , 

1 + 2x 2 > . (23) 



which are identical to the constraints derived in [14]. Note that the M2-brane with 
a single angular momentum is always stable, but that with two angular momenta or 
more an instability is possible. (For example, when one approaches criticality along 
a special directions of n equal non-zero angular momenta the critical values become 
x c ~ 2.118431 and x c ~ 2.193567 for n = 2 and n = 3, respectively.) 



B. Critical Exponents 



We will focus on the grand canonical ensemble. Critical exponents characterize the 
behavior near a critical point of the grand potential, also called the Gibbs potential, 
and defined as S = E — ST — J2i ^Wi- Conformal invariance implies S = ~E. Using 
(8) and (9), one derives 



pp 



VN a p (p + 1)p +1 

=-ct p+i h(i 



rpp+l'_ 



S P+1 

eP 



s T ) 



( 



p+i 



x !) 



3-D 



3-D + 2J2 



(24) 
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where C = (27T7) P+1 . One can also show that the inverse temperature of the 

solution is 



e = l 



8irp 



t (D-i)(p + iy 



( 



3-D 



3-D + 2J2 



(25) 



The expression (24) for £ is in terms of the x i: and it is straightforward to make a 
Taylor expansion around a point x ic on the boundary of stability. There is no singular 
behavior in this expansion. However, if £ is expanded around the chosen point x ic 
in powers of fli — f2 ci , there are in general fractional critical exponents, 7^. This is 
possible when in the expansion 



X, — LO ic = 27C I Aij 



u a d dy H 
where = - — — , 



X]c) + dx h 



\Xj ^jr'c) {% k %kc) ~\~ 



(26) 



the first term on the right hand side vanishes. There will always be directions of 
approach in Xj-space where this term does vanish: has vanishing determinant at 
the boundary of stability, and if we set Xi = Xi C — rjVi where AijVj = 0, then 



Ui - u ic ~ c 2 r] 2 + c 3 r] 3 + ... 



c ~ c 2 r] 2 + c 3 r] 3 + 



(27) 
(28) 



for some coefficients q and q. There is no linear term in 77 in (27) by construction of v «. 
There is also no linear term in (28) because such a term would cause some component 
of the angular momentum to become infinite as rj — > 0, and this we can rule out 
with calculations in the microcanonical ensemble. Eliminating rj one can obtain an 
expansion for £ of the form 



£ ~ (analytic) + \u>i 



1 2-7 



+ . . . 



(29) 



where the term shown explicitly is the leading non-analytic term. In the absence of 
numerical coincidences such as the rf term vanishing in (27), or the relative coefficients 
between r/ 2 and rf being identical in (27) and (28), we will have 7 = 1/2. It is important 
to keep in mind that the exponent 7 is a function both of the particular point on the 
boundary of stability one is approaching and of the direction of approach. In the cases 
we have examined, generically 7=1/2 with some isolated exceptions which we will 
now describe. 

Let us take set n angular momenta equal, so that x\ — . . . — x n — x, and send x 
to its critical value x c with the other Xi vanishing. In this case the expression for the 
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common value uo of uoi = in terms of x can be inverted explicitly, and it takes the 
form: 

* = l^TD ± [(77)' + (2n + 3 - £>)(£> - 3)] 1 / 2 ) , (30) 

Half integer powers appear in the expansion of x around Q — Q c when (2n + 3 — D) (D — 
3) < 0, and this leads to 7 = 1/2. This situation obtains when n — 1 for all D, and 
also for n = 2 in the case of the M2-brane. In fact n — 1 is essentially the generic 
case: approach to the boundary along a normal direction away from any special point. 
When n > 1 for the D3-brane or the M5-brane, and when n > 2 for the M2-brane, 
(30) leads to an expansion of £ which involves only integer powers of uj — oo c . 

On account of (26), there can be as many independent directions where 7 is fractional 
as there are independent vectors annihilated by A^. In order to determine when A^ in 
(26) has multiple zero eigenvalues, an efficient method is to determine the entire locus 
of points where the polynomials in (20) vanish. Let us call one of these polynomials 
p(xi), and note that p = Ai . . . X n up to a factor which is never 0, where Aj are the 
eigenvalues of Aij\ Xk=Xkc . The locus p = is the union of the vanishing sets of the 
Aj. At a generic point on p = 0, only one Aj = 0. Near a point of p = where 
the vanishing set of one eigenvalue crosses the vanishing set of another, the equation 

p = does not define a manifold; rather, it defines two or more manifolds crossing. 

I D-1 I 

We can find these points by solving p = and dp = J2i=i J^ 6 ^ = simultaneously. 
For the M5-brane, there is no solution: / = truly defines a manifold, and the two 
eigenvalues of A^ never simultaneously vanish. For the D3-brane and the M2-brane, 
the only solutions are where all x\ = 1. For the D3-brane, this is a branch point where 
two sheets of p = cross, and correspondingly two eigenvalues of Ay vanish. For the 
M2-brane, three sheets of / = cross, so three eigenvalues of Ay vanish. In general, 
the number of eigenvalues vanishing at a branch point is indicated by the lowest degree 
term in the Taylor expansion of p around the branch point: that is, the order of the 
branch point. 

C. Field theory analysis 

In [1] a field theory analysis of spinning D3-branes was attempted, using the fact 
that transverse angular momentum on the brane amounted to R-charge. In this section 
we will extend the analysis from one angular momentum to three. We will also briefly 
consider the M2-brane and the M5-brane. 

The on-shell degrees of freedom in a single D3-brane's world-volume theory are 
two vector boson states (counting positive and negative helicity separately), six scalar 
states, and eight fermion states. The R-charges of these states are specified by vectors 
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in the weight lattice of 5*0(6). The on-shell degrees of freedom for a single M2-brane 
are eight scalar states and eight fermion states, with R-charges specified by weight 
vectors of 5*0(8). The on-shell degrees of freedom for a single M5-brane are three 
self-dual two-form states, five scalar states, and eight fermion states, with R-charges 
specified by weight vectors of 50(5). All these states are tabulated in equation (31): 

{vector: 2(0,0,0) 
scalar: (±1, 0, 0) + (0, ±1, 0) + (0, 0, ±1) 
fermion: (±±,±±,±±J 



on-shell M2-brane states, 
50(8) quantum numbers 



scalar: (±1, 0, 0, 0) + (0, ±1, 0, 0) 

+ (0,0, ±1,0) + (0,0,0,±1) 

fermion: (±±, ±±, ±|, ±0 , 

even number of - ^-q 



{two-form: 3(0, 0) 
scalar- (±1 0) ± (0 ±1) ± (0 0) 

fermion: 2(±±,±±) 

For p = 2, 3, or 5, let a be a vector which runs over the sixteen p-brane weights 
given above, and set Sq = 1 when the corresponding particle is a boson and S3 — — 1 
when it is a fermion. We will calculate in the grand canonical ensemble with "volt- 
ages" Qi, f2 2 , Q[(d-i)/2\ corresponding to the elements of the Cartan subalge- 
bra of the R-symmetry group (for instance, in the case of D3-branes, we have Qi, 
Q2, and ^3 for the three elements of the Cartan subalgebra of 50(6)). Denoting 
uj = (u>i, u>2, ■ ■ ■ , w\(d-i)/2\ ) where uji = fli/T, we have in analogy to (26) of [1] an 
expression 

i = T I Wr^ Ss ' 0g(1 " s «- e """ 1+M ) = vCgya, I) L. !>+1 ( Sa e-) (32) 

for the grand potential per unit volume, £, which is also the pressure. As in the case 
of only one angular momentum, the integrals diverge for any real nonzero u because 
there are charged massless bosons. The integrals in (32) are all special cases of 



L 



dk k n ~ l log(l - e^ k ) = -r(n) Li n+1 (e^) , (33) 





which is carried out assuming complex \x and using analytic continuation. The difficulty 
in defining the integrals for real Q is apparent from the branch cuts across the positive 
real \i axis of Li p+ i(e M ). In [1] it was proposed to define the integrals using the principle 
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value prescription, which in this case amounted to replacing the polylogarithms by their 
real parts. This prescription together with the identity 



Ln/2J or(2i) m 



leads to 



z \ i / % 



(35) 



for the D3-brane, and to 



(36) 



for the M5-brane. For the M2-brane, the grand potential is transcendental, and cannot 
be simplified substantially beyond (32). At small uJi one has the expansion 

For N coincident D3-branes branes in the supergravity limit, the world- volume theory 
is strongly interacting. We nevertheless suggest (35) as a first approximation to the 
grand potential of 3 + 1-dimensional M = 4 SU (N) gauge theory at high temperature, 
where now we have scaled a factor of iV 2 out of £ as in section II A. The specula- 
tion is that even though the fundamental excitations are strongly interacting, for the 
purposes of thermodynamics we can think of 0(N 2 ) quasi-particle excitations which 
have the same multiplet structure as the fundamental degrees of freedom. For the M2- 
brane, agreement with supergravity requires 0(N 3 ^ 2 ) quasi-particles, whereas for the 
M5-brane, 0(N 3 ) are required. These peculiar scalings are borne out by absorption 
calculations [16], but are otherwise poorly understood. 
For the rest of this section we will focus on the D3-brane. 

Thermodynamic stability requires that £ be a subadditive function of T and the 
Because of the scaling laws due to conformal invariance, this amounts to being able 
to express the u>i as single- valued functions of the ratios ji/s. It is straightforward to 
check that this is possible on the region of cu space containing the origin where 

„ dg dg dg d 2 g dg d 2 g \ ^ n , ooX 

ouJkOLVf ouJk Otoe ouJk ouJiOuJe ouji ouokOuj£ I 
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In (38), which generalizes (42) of [1], the index % is summed over in the last two terms, 
and k and I are the indices of the 3x3 matrix whose determinant is taken. 

The mean field theory analysis of [1] also generalizes in a simple way. We will omit 
some details which do not depend on the number of angular momenta. The basic idea 
is that there is some "self-field" contribution to the voltage, which we express as 

if = u + f(j/T 3 ) . (39) 

Here tp is the total voltage, u is the applied voltage, and f{j/T 3 ) is the self-field contri- 
bution, given in terms of an unknown function /, and in terms of an R-charge density j 
obtained by integrating over the thermal occupation factors (using the principle value 
prescription) in the presence of the total voltage tp. For notational simplicity we have 
dropped the arrows from ip, u, and j, but they are indeed three-dimensional vectors, 
and derivatives with respect to them are also three-dimensional vectors. Thus the per- 
meability n = dtp/duj is a 3 x 3 matrix which can depend on uo. The simplest extension 
of [1] is to assume f(x) = —x, which implies that the permeability at zero voltage is 
/i — 3I. That is, the D3-brane responds diamagnetically to an applied voltage. We 
will show how this ansatz reproduces the critical behavior found in supergravity. 
Mean field theory leads to a parametric expression for the free energy: 

-X = aH =g(<p) 

d a < 40 > 

LU = (P+ — . 

dip 

The first line in (40) uses (32) with u replaced by tp. The second line is a rephrasing of 
(39) using the ansatz f(x) = —x. The function g(tp) is the same as appeared in (35). 
Since g{tp) is analytic, the only way nonanalyticity can arise is through inverting the 
relation between u and tp. Nonanalyticity can arise in this context only if the inverse 
permeability, 

i du „ d 2 g 

" =aT 1 + a?' (41) 

acquires a zero eigenvalue. Thus the boundary of stability is determined by det = 0, 
which is equivalent to 

576tt 6 - 48tt 4 (v?i + tfj + tp\) - 207r 2 (y^ + tp\ + tp\) 

+ ^\tp\tp\ + tp\ip\ + tpltpl) - (tp\ + tp 6 2 + tp\) + Gtpftpltpl (42) 
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This complicated expression reduces to \ifi\ < 2tt when if 2 = y?3 = 0. It exhibits the 
same qualitative behavior as (20), in that the boundary of stability is a generically 
a single zero of the left hand side, with the sole exception of all ipj equal, when the 
boundary of stability becomes a double zero. As before, this indicates that generically 
only one eigenvalue crosses through at the boundary, but for all angular momenta 
equal two eigenvalues vanish at the same time. We obtain the same critical behavior 
from the mean field theory analysis as from the supergravity, and for basically the same 
reasons: making the expansions 

u - u c = (l + g"(if c ))(if - f c ) + \g'"{f c ){v - y? c ) 2 + • • • 

t ( 43 ) 

= g(f c ) + g'(<Pc)(<p - Vc) + y"((f c )((f - fc? + ■ ■ ■ , 

we see that 7 = 1/2 for variations of u in the directions in which has a zero 
eigenvalue, while 7 = in the other directions. The algebra is the same as in (27) 
and (28), except that there 77 represented a deviation parametrized in terms of x — x c , 
whereas in (43) it would be a deviation parametrized in terms of ip — ip c . 



III. PHASE MIXING 



So far our goal has been to demonstrate that spinning branes with a sufficiently high 
angular momentum become thermodynamically unstable. We have not answered the 
obvious question: what form does the instability take? In this section we will entertain 
two possibilities: 1) the branes fragment and fly apart in the transverse dimensions; and 
2) the angular momentum density localizes on the branes. For the sake of specificity, 
we will focus exclusively on N coincident D3-branes with only one angular momentum. 
In the notation of [1], 

27tt 2 {j/vy 1 

X 8N 2 {E/Vf > 3 1 ' 

puts us outside the region of stability as calculated in the grand canonical ensemble. 

The naive expectation is that D3-branes with x > 1/3 vigorously shed their angular 
momentum by radiating either closed strings or individual D3-branes. In an examina- 
tion of the radial potentials for closed string emission and for equatorial motion of a 
test D3-brane in the background of N spinning D3-branes, we found no evidence for 
a special unstable mode developing at x — 1/3- This is not necessarily a problem: as 
commented upon in the introduction, the nature of the instability is not so much that 
there is some runaway mode of the system by itself, but rather that it is impossible to 
establish equilibrium between it and a heat bath. 
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The following simple estimate suggests that for \ > 1/8 D-brane fragmentation is 
thermodynamically allowed. 3 Suppose the the N D3-branes split into k equal frag- 
ments. As the fragments fly apart, they can acquire orbital angular momentum. As 
an idealization of this complicated dynamical fission process, let us suppose that most 
of the angular momentum winds up in the orbital motion of the fragments, and that 
most of the energy winds up in their non-extremality. 4 As a further idealization, we 
will assume that the process occurs uniformly across the D3-brane world-volume. By 
compactifying the world- volume on a small T 3 one could ensure such uniformity. From 
[1] we take the formula 

2 5 /4 3 -3/4 v /^ A 7 e 3/4 

s(e,j;N)= (45) 

for the entropy per unit world- volume in terms of the energy and angular momentum 
per unit world-volume. Since we will be considering different values of N, we have 
reverted to the conventions of [1]: e = E/V, j = J/V, and s = S/V, without the 
factors of iV 2 in the denominator which we used in section II A. The ratio of the 
entropy of the final state (k identical non-spinning fragments) to the initial state (a 
single cluster with spin) is 



/TT^ + Vx. (46) 



8f_ = ks (f>°; f) _ 

Si s(e,j;N) ^ y/k 

The process is entropically favored if Sf > Sj, which is equivalent to \fk < y/1 + x+ \fx- 
Kinematics forbids the k = 1 process, so we get the lowest critical x from k = 2 (two 
fragments): x > 1/8, as advertised. With k — 3 (three fragments) we get x > 1/3, 
which coincident ally corresponds exactly to the boundary of the domain where s(e,j) 
is subadditive. Making the fragments unequal only increases the bound. For very large 
X, large values of k will be accessible, and we interpret the fragmentation as a Jeans 
instability of the disk of D3-branes described in [9]. 

The above estimate does not convince us that brane fragmentation is the dominant 
instability for x slightly larger than 1/3. We would be more optimistic if an instability 
of the horizon could be identified. In any case, brane fragmentation would violate 



3 We thank F. Larsen, E. Martinec, and P. Kraus for a stimulating discussion in which this 
calculation was suggested to us. 

4 If we accounted for the energy that goes into Hawking radiation and the motion of the 
fragments, the bound would be x > Xo for some xo > 1/8. 
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singularity theorems, so we may expect it to have a very long time-scale in the classical 
limit of many coincident branes. 

We now turn to the second possibility: that for x > 1/3 the angular momentum 
density j tends to localize on the world-volume. If they are thermodynamically per- 
mitted, inhomogeneities in j might nucleate locally on a shorter time-scale than brane 
fragmentation. In fact, the time scale should be on the order of the inverse temperature 
since there is no other scale in the problem. 

Let us regard (45) as an approximate expression for the entropy density in terms 
of energy density e and angular momentum density j which are allowed to vary over 
the world-volume. Diffusion and interactions will tend to keep the temperature T and 
the "voltage" Q (properly the angular velocity at the black hole horizon) spatially 
constant. However, for a given (T,Q), there are two possible values for energy density 
and angular momentum density. (That there are two and only two possible values can 
be seen from the fact that the equation of state, (20) of [1], is do uble- valued) . One of 
them, call it (e, j), will lie in the stable region where s is subadditive; the other, call 
it (e, j), will lie in the unstable region. Thus we are led to consider mixed states of 
the two phases (e, j) and (e, j). We will argue that such states are entropically favored 
over states with uniform (e, j) for x ~ 0.325, although only for x > 1/3 can a uniform 
state evolve smoothly into such a mixed state. The mixed states are only saddle points 
of the entropy functions, not maxima, because they involve an unstable component. 

We emphasize that the following calculation is not on the same footing as the usual 
analysis of phase coexistence, because one "phase" in the current context is unstable. 
Rather, it is intended as a preliminary probe of how inhomogeneous e and j may arise 
near x — 1/3, based on the idea that T and Q should remain at least approximately 
uniform while these inhomogeneities are first forming. Because we are relying on the 
known supergravity solutions as our guide to the thermodynamics, we are not equipped 
to address the questions of coherence length, surface tension, and stability of phase 
boundaries. Ideally, one would like to establish the existence of a second stable phase, 
perhaps one where the branes are separated in the transverse dimensions. Then a more 
standard equilibrium analysis would become possible. 

Assume that a fraction A of the spatial volume is filled by the unstable phase (e, j), 
and the remaining fraction 1 — A is filled by the stable phase, (e, j). Then the average 
energy density, angular momentum density, and entropy density are 



e, 



Ae+ (1 - A)e 



Jav 



Aj + (1 - \)j 



(47) 



s, 



av 



Xs(e,j-N) + (l-X)s(e,j;N) , 



18 



where s(e,j;N) is the function (45). With e av and j av fixed we want to extremize 
s a v by varying A. That is, in an isolated system where the total energy and angular 
momentum are fixed, we extremize the total entropy by distributing energy and angular 
momentum optimally between the two phases. This problem is mathematically well- 
posed because e and j can be determined as functions of e and j, and then the first 
two lines of (47) amount to two (non-linear) constraints on three variables, leaving A 
as the only free variable. 

More explicitly, for given (T, f2), there are two solutions to the equations 




1 = v ^F3 1 / 4 /2 3 / 4 1 j 

and if one of them is (e, j), then the other, (e, j), can be expressed as 

4 _ 




, (49) 

where w = ^ — — , X = X + -777^ — 7 • 

1 + X w(2 + w) 

With the help of (49) one can show that 

s(e,j;N)=s(e,j;N) {^j'' ^ ■ (50) 
Because of the overall scale invariance of the theory, it is equivalent to extremize 

(Jav = "-f (51) 

with fixed 

27n 2 j 4 

Xav = &W~ef v ■ (52) 

Using (47), (49), and (50), one can obtain a av and Xav explicitly as functions of x an d 
A. The domain of x is [0, 1/3] because it is by assumption the ratio pertaining to the 
stable component, and the domain of A is [0, 1]. We have verified numerically 5 that the 



'Plots and numerics were obtained using Mathematica. 
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contours of constant Xav can be written as single- valued functions, A = A Xov (x). Thus 
the entropy along a given contour is extremized when 

fd^A = f^\_fdg\ (*=^ = . (53) 

In figure 1 we show contours of constant Xcre cut by a dashed line indicating the solution 
set of (53). 

From figure 1 one can conclude that for Xav > 1/3 the entropically preferred state, 
among states with uniform temperature and voltage, is a mix of the stable and unstable 
phases, with A <: 0.46 and A — > 1 as Xav —> °o. The Xav ~ 0.322 mixtures are always 
entropically disfavored: for such Xav the preferred state has A = 0. There is a narrow 
band of values, 0.322 ^ Xav < 1/3, where there are two far-separated competing states 
(see figure 2): one has A = and one has A m 0.46. At Xav ~ 0.325 there is a 
"transition" where A jumps from to 0.46. The temperature and voltage also jump at 
this transition. 

We regard this as preliminary evidence that there is a first order phase transition at a 
value of Xav strictly less than 1/3. The critical behavior associated with the boundary 
of stability would be cut off finitely below the critical angular momentum by phase 
separation. A more definitive treatment would be possible if one could identify a ther- 
modynamically stable phase into which the system falls when the angular momentum 
density is sufficiently high. Our treatment has assumed that this new stable phase is 
very slow to form compared to the time-scale on which small inhomogeneities in j can 
nucleate. This is reasonable if the new stable phase is some version of a multi-center 
solution, because passing to such a phase would violate classical singularity theorems. 

In [14] it was proposed that in the canonical ensemble, stability persisted up to 
X ~ 1-26. As a general rule we agree that thermodynamic stability of a system depends 
on its environment. However the results of this section indicate that the instability of 
spinning D3-branes at x = 1/3 is not easily avoided by changing ensembles. Consider 
a state with total e and j arranged so that Xav is just slightly larger than 1/3. The 
pure phase state would correspond to where one of the grey lines in figure 1 intersects 
the top of the frame. At this point, x — Xav This is the sort of state that the authors 
of [14] regard as stable in the canonical ensemble (though not in the grand canonical 
ensemble). However, the system may increase its entropy continuously by flowing down 
along one of the grey lines, until it reaches the dashed black line, where entropy is a 
maximum at least among two-component mixed phase states of the prescribed Xav In 
other words, the putatively stable state represents only a saddle point of the entropy, 
not even a local maximum. This is in contrast with states with 0.325 ^ Xav < 1/3, 
where a system in a pure phase state can increase its entropy only by making a big 
leap, all the way along one of the fine black lines from the bottom of the frame to the 



20 



FIGURES 



1 


/////// i 




0.8 






0.6 






X 

. 4 










0.2 














0.05 0.1 0.15 0.2 0.25 0.3 



X 



FIG. 1. The solid black lines are contours of constant Xav with \ av < 1/3. The grey lines 
are contours of constant \av with \av > 1/3- The dashed line is the solution set of (53). 
The region below the dashed line is where entropy along a contour increases as one moves to 
larger x- If a contour intersects the dashed line, then the maximum of the entropy function 
along that contour is at the first intersection with the dashed line; otherwise the maximum 
is at the lower border A = 0. 




FIG. 2. The dimensionless ratio a av = s^ v /e^ v is plotted against x f° r two values of \av 
on either side of the phase transition: in a) Xav = 0.324 while in b) Xav = 0.326. In a) we 
see that the x = 0.324, A = 1 solid phase wins out entropically; in b) the mixed phase with 
X ~ 0.075 and A « 0.46 is entropically favored. 
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leftmost intersection of the fine black line with the dashed black line. In a nutshell, the 
pure phase states represented by the bottom of the frame in figure 1 are metastable 
against separation into a two-component mixture, but the pure phase states represented 
by the top of the frame are not even metastable. These statements are independent of 
the environment because we are talking only about shifting energy around within the 
world- volume of the D3-branes. The key point is that D3-branes have spatial extent, 
so any little bit of it can be regarded as a "system" in thermal contact with a "heat 
bath" consisting of the rest of the brane. That is what makes the grand canonical 
ensemble relevant. With black holes the story is different: because there is no spatial 
extent, it is doubtful that one can regard one part of the black hole as a heat bath for 
another part. If the D3-brane world- volume were compactified on a scale smaller than 
the thickness of a domain wall between phases, then the situation would be physically 
indistinguishable from a black hole, and calculations in the canonical ensemble might 
become relevant. 

IV. EUCLIDEAN SPINNING BRANES 
A. Spin and thermal boundary conditions 

Euclidean black brane solutions can be obtained from the Minkowskian ones by 
sending t — > —it and £j — > i£j, where the £j are the angular momentum parameters 
that enter into the solutions of [8]. Complexifying them is necessary in order to keep 
the dtd(f) parts of the metric real. The consequences for the thermodynamics can all be 
reasoned out from the fact that the angular momentum becomes complex. In brief, one 
can retain equations (3) through (14) with the replacements i/j — > iyj and Xj — > ixj, 
but uh — * Uh with no factor of i. These replacements have a profound effect on the 
physics: as we shall now argue, the phase space is \xj\ < 1, and points on the boundary 
correspond to supersymmetric periodicity conditions in the Euclidean path integral. 

The new version of the first equation in (14) is 

Vh,E = 11(1 - *?) • (54) 
j 

This relation becomes singular when any of the Xi reaches 1 in absolute value, so the 
phase space is indeed the cube \xi\ < 1. The new version of (11) and (13) is 

-f = 2m^ = -2ni — J — ^ . 55 

% 1 - X te-D-Xk£% 

We now want to study the limits of (55) as we approach the boundary of the cube. 
Consider fixing some number of the Xj at values finitely far from the horizon and then 
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letting all the others approach the horizon at the same rate: for these others, say 
Xj = Sj(l — e) where Sj = ±1 and e — > 0. For convenience define Sj = for all the Xj 
whose values are fixed in the interior. Let the number of Xj's which are approaching the 
boundary be n. Approaching a generic point on the boundary corresponds to n — 1; 
n > 1 means one is approaching a point on an edge or a corner of the cube rather than 
a face. It is simple to see from (55) that 

flj 2m . . 

T n 3 y ' 

For example, if we send x± — > 1 and hold all the other X{ fixed, then Q±/T — > 27ri and 
all other fy/T -> 0. 

To understand what is happening in the world-volume theory, let us first recall that 
in the absence of rotation, the usual thermal partition function, Z = tre _/3i? , can be 
computed by a Euclidean path integral with period 1/T for Euclidean time. Boundary 
conditions on fermions are anti-periodic. To add real angular momentum — the kind 
that we have studied in all the previous sections — one would introduce a chemical 

— * 

potential f2 for particles with R-charge: acting on a state \a) whose R-charge is specified 
by the vector a in the weight space of the R-symmetry group, the hamiltonian H would 
be replaced by H—a-Q (cf. (32)). But in this section we are studying Euclidean angular 
momentum, which corresponds to pure imaginary Q. Thus 



e 



^/T\&) = e™' T \a) = R\a) , (57) 



where H = (H 1 , H 2 , . . . , H^ D _ 1 y 2 \) is a basis of generators for the Cartan subalgebra. 
R = e H ' n / T is an element of the compact R-symmetry group SO(D — 1) (or more 
precisely its covering group) when Q is pure imaginary. The net effect of turning on Q, 
as we see from (57), is to insert R into the partition function: instead of Z — tre~^ H , 
now we have Z = tr Re~^ H . In the Euclidean path integral, R just specifies a twist in 
SO(D — 1) which one performs on all the fields before identifying from ts = 1/T to 
ts — 0. This twist modifies the usual (thermal) boundary conditions according to the 
R-charge of each field. 

If we set njT = 2m and all other Qj/T = 0, then R = (-1) F , where F = for 
bosons and 1 for fermions. That is, we recover supersymmetric boundary conditions. 
Maximal Poincare supersymmetry is recovered in this case because Rrj = —t] for any 
spinor rj of SO(D — 1). This does not mean that conformal invariance is recovered. In 
fact it is not: the near- horizon geometry is not anti-de Sitter for Euclidean spinning 
branes, even in the large spin limit we are considering. We are led to conclude that 
the geometry must be describing a physical state analogous to those considered in [9], 
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where Higgs VEV's break conformal invariance. 6 

For n > 1, we can judge whether any part of supersymmetry is restored by deter- 
mining whether there are any spinors rj of SO(D — 1) which satisfy Rrj = —r], where 
R is the element of the covering group of SO(D — 1) specified by Qj/T = l-nijn for 
1 < j < n- It is straightforward to verify in this way that a fraction (l/2) ra_1 of 
supersymmetry is preserved in all cases except n = 4 for the M2-brane, where 1/4 of 
supersymmetry is preserved. 

By studying Killing spinor equations in the spinning brane supergravity geometries 
one can confirm the presence of unbroken supersymmetry. This was essentially done 
in [18] for the M2-brane. To be more precise, [18] includes an investigation of the 
Killing spinor equations for BPS-saturated R-charged black holes of N = 8 D = 4 
gauged supergravity. The results of [19,13] indicate that Kaluza-Klein reduction of 
spinning M2-branes leads precisely to large R-charged black holes of the type studied 
by the authors of [18]. They found that 2 5 ~ n supersymmetries are preserved by black 
holes with n — 1, 2, or 3 charges, while adding a fourth charge does not break any 
additional supersymmetry. This counting is in agreement with the previous paragraph: 
for example n = 1 corresponds to 16 supersymmetries, which is the maximum when 
conformal symmetry is broken by the physical state. 

The weight vectors listed in (31) can be used to determine the allowed momenta fog 
around the Euclidean S 1 . These momenta determine tree- level masses in the dimen- 
sionally reduced theory. For n < 4 there is Bose-Fermi degeneracy in all cases (M2, 
D3, and M5). (See also figure 3.) The M2-brane with n = 4 exhibits no Bose-Fermi 
degeneracy. This may be an indication that our description of the world-volume theory 
of many M2-branes in terms of excitations with the same SO (8) quantum numbers as 
for a single M2-brane is too simplistic to capture even the rough outlines of the field 
theory. 

In order to further shed light on the nature of supersymmetry restoration in field 
theory and to make comparison with the supergravity results, let us turn to the study 
of the thermodynamic quantities. 

We obtain flj/T = (2iri)/n exactly when Xi = ±1 (i — 1, • • • , n), but other thermo- 
dynamic relations become singular in this limit. So we choose uii = (±)27ri/[n(l — e)] 
(xi = ±(1 — e)) (i = 1, ...,n, e 1), and for the sake of simplicity we set other u/s 
to zero. This is still not the most general way one could imagine approaching the 
boundary of phase space: for instance, one could set Xi = ±(1 — r^e) with arbitrary but 
fixed r\. This does not lead to any more general scaling behavior than what we will 



6 Related ideas have recently been explored in [17], although there are subtleties with the 
Wick rotation in comparing with the present work. 
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observe, but the coefficients in equations (58) through (67) and (71) would change. 

First we compute the grand potential for the field theory obtained from (32) with 
the super symmetric particle content given in (31). The field theory results in the e — > 
limit are as follows: 

D3: ^ = -7T 2 (§e 3 , |e 2 ,±e ), n= (1,2,3), (58) 
M5: JL = -vr 3 (±e 4 , ^e 2 ) , n=(l,2) , (59) 
M2: = - (7rln2 e 2 ,vr(| + hi^)e 2 ,2.0988e , 0.33480 ) , n= (1,2,3,4). (60) 

For a single M5-brane, or D3-brane this field theory analysis is sound: it amounts 
simply to counting states in a free field theory. (It is not even plagued by the divergences 
that we used the principle value prescription to solve in section II C: those divergences 
were for uii real, and are completely avoided by pure imaginary For N > 1 our 
understanding of the field theory at finite temperature is limited. The naive approach is 
to hope that it has some properties in common with 0(N 3 ) copies [0(N 2 ) copies] of the 
N — 1 theory in the M5- [D3-] brane. This approach is not well motivated physically, 
but it does give the correct scalings with iV and e as predicted by supergravity. The 
field theory on multiple M2-branes not very well understood. It is based on a particular 
infrared limit of the D2-brane gauge theory where 5*0(8) symmetry is thought to be 
recovered (see for example [20,21]). For a single brane, a simple Hodge dualization of 
a vector field into a scalar field suffices to show this, but for multiple branes there does 
not seem to be a simple argument directly from field theory. 

The positive powers of e in (58)-(59) appear because of the Bose-Fermi degeneracy in 
the spectrum of allowed Ice- This is substantially the only indicator of supersymmetry 
in this free field theory computation. 

Let us now compare the field theory results to those obtained in supergravity. There 
the grand potential £ is given by (24). In the e — > limit, £ has the following scaling 
behavior: 



£ „ TP+h^ +p+1 . (61) 



The specific results are: 



D3: = -7T 2 (l6e 3 , |e 2 , ) , n= (1,2,3), (62) 



M5: ^ = -vr 3 (f e 4 , ^e 2 ) , n = (l,2), (63) 



^6 



M2 : jr 3 = -- 2 (f , ^e 2 , §f e 3 / 2 , ) , n = (1, 2, 3, 4) . (64) 
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It is instructive to trace the behavior the inverse temperature in this limit. Using 
the expression (25) one obtains: 

D3 : (3 

M5 : 13 
M2 : 

Note that for finite temperature to remain finite r ~ (e 1- ? : ^Hr- ,e^) for D3-, 
M5- and M2-brane, respectively. 

Comparing the supergravity results (62) through (64) with those the field theory 
analysis ((58) through (60)) reveals that except for rational prefactors the supergravity 
and field theory results are in good qualitative agreement for both D3- and M5-brane, as 
expected since a single D3- and M5-brane field theory analysis is believed to be sound. 
On the other hand there is a qualitative discrepancy between the two approaches for the 
M2-brane case, again indicating that the field theory of M2-brane is poorly understood. 

Most importantly, we see that the energy density e = — p£ scales in all cases with 
a positive power of e. Recall that e is the energy density over and above the energy 
density of a BPS-saturated brane configuration. By the same argument that zero 
energy states in a supersymmetric theory are supersymmetric states, we conclude that 
a state with e = preserves at least some of the original supersymmetries of the brane. 

We would like to conclude this subsection by commenting on the supergravity anal- 
ysis of the spectrum of the theory. This will become relevant in our discussion in the 
next subsection of models of confinement based on spinning branes. Let us concentrate 
on a minimally coupled scalar 0, that is a scalar whose linearized equation of motion 
is the Laplace equation in the background of the Euclidean spinning brane: 

-^=d,{^Tgg^)d^ = o. (68) 

With the ansatz = e lkx x(f), where x is a world- volume direction, and (for the sake of 
simplicity) x( r ) depends only on the radial coordinate r, one obtains a radial equation 
of the following form: 7 



7 In deriving the above equation the full metric of the spinning brane has to be used. We 
used the metric of the general D-dimensional rotating two-charge black hole solutions given 
in [8]. Again, the corresponding spinning (D3, M5, M2)-branes are obtained by turning off 
one of the two charges of the (D = 7, 6, 9-dimensional) black hole, and lifting the solution to 
ten and eleven dimensions, respectively. 




(1,2,3), (65) 

n= (1,2), (66) 
(1,2,3,4). (67) 
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D odd 



D even : 



rR D-3 
rR D-3°r 



-1 



LVJ 

] J (r 2 - £ 2 ) - 2mr 



i=i 



9 r x(r) = k 2 x{r) , 



n ( r2 - ^) - 2mr 

i=i 



^x(r) = A; 2 x(r) . 



(69) 



These wave equations were derived in [10] and [11] for the D3-brane and the M5-brane, 
respectively. In the field theory, the values of \k\ will be masses of ++ glueballs. 

In terms of the notation introduced in this paper these wave equations assume the 
form: 

LVJ 

//;; ;! -<- 1 n (* 2 -^)-* 



D odd : x 1 d T 



i=l 



dxx = y H ( kr oY 



R 



r 



D-3 



X , 



D even : x 1 d T 



i=i 



d xX = y 2 H (kr y ^ \. (70) 



Our radial variable is x = The quantities Xi, i/i, yn, ?"o, and R were introduced at 
the beginning of section II A, but because we are now operating in Euclidean space we 
must recall the replacements and yj — >■ iyj. When n equal Euclidean angular 

momenta becoming large at fixed energy (that is, Xi — > ±(1 — e) for i = 1 . . . , n), the 
limiting behavior for k is 

^M glucball = ^Q(2e) 1 -tT, (n = 1, • • • , ^) , (71) 

where T is the Hawking temperature (25) (which at the boundary takes the values 
(65)-(67)), and Q is a pure number determined by the following eigenvalue equation: 

LVJ 

x- 1 n (* 2 -^ 2 ) 



D odd : x 1 d r 



i=i 



D even : x 1 d r 



n ^--f) 



dxX = -Q X, 



d xX = -Q 2 x, 



(72) 



where now x 2 = 1 (i — 1, . . . ,n) and other x/s are set to zero. (As before, with non- 
zero x 2 7^ l's and Xi = ±(1 — r^e), the equations (72) are of the same form, and the 
mass gap (71) has the same scaling behavior with e, but the actual numerical prefactors 
in (71) change.) Discrete eigenvalues Q arise from demanding that x is normalizable. 8 



Certain cases of (72) can be solved in terms of hypergeometric functions — for instance 
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Thus Mg a p ~ -JIT for n = 1, M 9fip ~ T for n = 2, and M gap diverges for n > 2. For 
n > 2 it is possible that there is actually no discrete spectrum in the e — > limit. This 
point is under investigation. 

The mass gap can also be quantified in terms of the ratios of the angular momentum 
densities at the boundary, ji, (i — 1, • • • ,n) and e. It is straightforward to obtain the 

n(p+l) , 

following scaling behavior: ji ~ T p e ' A - D . One can easily show that 

f 1 2 \ 1/(p - 1) , n( P+ l) , „ 

M gap ~ ( =j J ~ Te 1- ^-^- 1 ) = Te 1 " 2 . (73) 

Some of the observations we have made about the scaling of the mass gap and 
about supersymmetry restoration also follow easily from the results of [22,6,10,11]. 
The connection that does not seem to have been made there is that the field theory 
itself recovers supersymmetry in the limit of small e because of the insertion of an 
element of the R-symmetry group into the partition function. 

B. Models of confinement 

The usual route to QCD 4 [23] is to compactify the M5-brane worldvolume on a 
circle twice. The first time, one uses supersymmetric boundary conditions on a circle 
of circumference 2nRi, and by the M-IIA duality the result is a D4-brane which, like 
the M5-brane, is near-extreme. The second compactification is applied to the Euclidean 
time coordinate of the D4-brane: it has period 1/T, and we impose thermal rather than 
supersymmetric boundary conditions. This second compactification is the important 
step. The five-dimensional theory one is starting with on the D4-brane is the Euclidean 
version of the SU(N) gauge theory with maximal supersymmetry which descends from 
the ten-dimensional M = 1 supersymmetric Yang-Mills theory. The resulting theory 
in four dimensions will also have SU(N) gauge fields, but no supersymmetry, and it 
is the candidate for QCD 4 . Tracing through the M-IIA relationship determining the 
string coupling in terms of the compactification radius R±, and then the relationships 
between the string coupling, the D4-brane gauge coupling, and the four-dimensional 
QCD 4 coupling #y M , one finds g\ M = 2^R{T. The relevant points for our analysis 
are: 1) the compactification process is unchanged when there is angular momentum in 
the directions orthogonal to the M5-brane; 2) the 't Hooft coupling A = gy M iV needs 
to be large for the supergravity to be valid. 



n = 1, D = 7 — and exact expressions for the eigenvalues Q were obtained in [17]. Some other 
cases can be solved in terms of Heun functions. 
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When there is no angular momentum, the fermions are anti-periodic in the Euclidean 
time used in the D4 — > QCD 4 compactification, so in the four-dimensional lagrangian 
they have masses ttT. It is usually assumed that at large 't Hooft coupling, the scalars 
pick up comparable masses. From supergravity one learns that the mass gap is also of 
order T. This is sensible from a field theory point of view: in the renormalization group 
flow of a gauge theory in which the 't Hooft coupling is large at energies comparable to 
the masses of the matter fields, one expects confinement at only lower energies. This 
is to be contrasted with a gauge theory where the 't Hooft coupling is weak at energies 
comparable to the masses of the lightest matter fields. There we have the standard 

one-loop relation M gap ~ exp (y—^^j ^matter, and b = yiV. 9 

When there is a large angular momentum, or, in the parlance of section IV A, small 
e, the story is rather different. The fermion masses are not all on the order of the 
temperature: as indicated in figure 3, some of them are on the order of eT. 
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FIG. 3. The spectrum of momenta in the compactified S 1 direction of the D4-brane world 
volume. The larger ticks indicate bosons; the smaller ticks indicate fermions. The first line 
follows from standard thermal boundary conditions. The second and third lines show what 
happens when one approaches the boundary of the phase space (large angular momentum for 
fixed energy), either at an edge (n = I) or a corner (n = 2). The splittings in these cases are 
on the order eT. In the close groupings of ticks in the second and third lines, the degeneracies 
follow the ratios 1:4:6:4:1 and 1:2:1, respectively. 

Supersymmetry is restored in the limit e — > 0. In the n = 1 case, for energies much 
smaller than T, only the center grouping of modes is available, and the theory is four- 
dimensional M = 4 super- Yang-Mills. In the n = 2 case, the center group of states fills 
out the Af = 2 super- Yang-Mills multiplet. 



3 We thank O. Aharony for pointing out to us the comparison with the one- loop analysis. 
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Supergravity calculations which we summarized in the previous subsection indicate 
that the scaling of the mass gap is M gap ~ ^feT for n — 1 and M gap ~ T for n — 2. 
Near the end of this section we explain why this scaling is rather at odds with our 
field theory expectations. For now we will regard the scaling of M gap as a non-trivial 
prediction of supergravity. 

In view of the existence of matter fields comparable to or much lighter than the 
scale of confinement, it is no surprise that large angular momenta did not suffice to 
decouple the "glueballs" charged under the global symmetry group which rotates these 
quarks among themselves [6] . The existence of such "Kaluza-Klein cousins" has been a 
persistent feature of all supergravity models of QCD 4 , but it need not be regarded as an 
artifact of the construction which must disappear in a full string theory treatment. In 
any confining gauge theory (including real-world QCD) in which there are fundamental 
matter fields with some flavor symmetry and with masses comparable to or lighter than 
the confinement scale, one expects to have flavored hadrons as well as flavor-neutral 
glueballs. It may not be reasonable to hope that a full string theory treatment of these 
backgrounds will decouple all matter fields and leave us with pure glue. It seems more 
likely that it will leave us with confining glue coupled to light matter. 

It makes sense that glueball mass ratios should be roughly the same [11] when one 
makes one angular momentum large as when one makes two large in any ratio except 
unity. In terms of the x« variables discussed in the section IV A and defined in (12), 
an M5-brane with one angular momentum large with the other zero corresponds to 
xi — > 1 while X2 = 0, whereas two large angular momenta with fixed ratio corresponds 
to xi — > 1 with X2/X1 = rj. As long as r\ 7^ 1, one is approaching the edge rather than 
the corner of the square which represents the phase space of the spinning brane, and 
the pattern of tree-level masses shown in figure 3 obtains. Glueball mass ratios should 
be and are approximately independent of r\ because the matter fields whose tree-level 
masses are small in comparison with the confinement scale. When 77 = 1, the pattern 
of tree-level masses is as shown in the n = 2 line of figure 3. We expect that glueball 
mass ratios will change at this point because the relevant matter fields have changed 
their tree-level masses. 10 

One difficulty with the original proposal of [23] is that the scale of the fifth dimension 
is the same as the scale of confinement. Thus any process whose energy is high enough 
to probe the gauge dynamics non-trivially — in the sense that colored states, not just 



The term "glueball" is somewhat of a misnomer since we are thinking of states which should 
involve adjoint matter fields even if they are singlets under the global flavor symmetry that 
the matter fields carry. The term "glueball" is in use because of the expectation (unlikely in 
our view) that the supergravity models discussed will correspond to pure QCD 4 . 
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an effective theory like soft pions, would be needed to analyze the process — would also 
probe the fifth dimension. One might therefore ask in what sense one has defined a 
four-dimensional gauge theory at any scale. The spinning branes construction is better 
in this regard because for n — 1 there is a large range of energies, y/eT <C E <C T, 
where the theory is four- dimensional and unconfined. The decoupling [6] of "glueballs" 
with Kaluza-Klein momentum in the fifth dimension is the tangible evidence that the 
length of the fifth dimension is much smaller than the inverse QCD scale. 

Other approaches to confinement avoid a fifth dimension altogether by using only 
D3-branes, either in type theory [24] or with a varying dilaton in type IIB theory 
[25] (see also [26]. These models have their own difficulties: in the former the gravity 
description breaks down in the ultraviolet, and there is a tachyon whose mass is on the 
string scale; while in the latter one seems to need some excited string state to start the 
dilaton flowing, and asymptotic freedom cannot be achieved. Both approaches involve 
naked singularities in the bulk. 

We summarize the taxonomy of holographic approaches to confining gauge theories 
in figure 4. The basic features of confinement are not expected to depend on the 
existence of matter fields or on their masses, provided that the beta function remains 
negative. All the gravity models have a string tension which is much larger than the 
square of the mass gap. This leads to some behaviors very unlike real QCD [27], in 
particular nearly flat Regge trajectories. In figure 4, we have indicated only the mass 
gap — that is, the mass scale seen by closed string probes. In all cases except 3), the 
QCD string tension would be located between this mass scale and the next highest 
scale; in 3) the relation of the QCD string tension and the temperature can only be 
determined given the relative sizes of e and the 't Hooft coupling g\ M N . 

Since only one large angular momentum is needed to make the scalings we have dis- 
cussed, it should be possible to construct variants of the spinning M5-brane approach 
by putting the branes at some spacetime singularity. Orbifold backgrounds would be 
expected to lead to bifundamental quarks rather than adjoints. There might even be 
chiral models analogous to those discussed in [28]. A model with bona-fide light fun- 
damental quarks could be difficult to construct, since the fundamental representation 
is the charge carried by one end of a string ending on a D-brane, and not typically by 
an open string with both its ends on the brane. 

There is one outstanding difficulty: why for n — 1 is the confinement scale much 
higher than the tree- level masses of the light matter fields (y^T ^> eT)? The most 
obvious resolution would be for at least some of the light masses to receive radiative 
corrections on the order yfeT from diagrams where heavy fields run around loops. 
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FIG. 4. Energies scales in various attempts to find a supergravity dual to confining gauge 
theories. Also included for purposes of comparison is a cartoon sketch of real-world QCD 
and a generic lattice model. The quark masses shown are the parameters that enter into the 
lagrangian: constituent masses. 



If this were to happen, then superconformal invariance would be broken at a scale 
on the order of y^T, and confinement would follow at an energy scale only slightly 
lower. But the following simple estimate suggests that perturbative effects cannot 
generate a sufficiently large radiative correction. Scalar mass renormalizations at one 
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loop typically arise from integrals of the form 



F(M) = J 



d 4 q 1 M 2 




(74) 



(2tt) 4 q 2 + M 2 8tt 2 



where we have used dimensional regularization and the minimal subtraction scheme. 
M is the mass of the heavy fields running around the loop. Typically, in a theory 
with supersymmetry broken by mass splittings for the heavy fields on the order AM, 
we have 5m 2 ~ (AM)F'(M) ~ (AM)Mlog^ on account of the near-cancellation of 
bosonic and fermionic contributions. With mass breakings such as the ones indicated 
in figure 3, the cancellation is even better: 5m 2 ~ (AM) M F iM \M) where Af is the 
degree of supersymmetry restored when AM = 0. In particular, 5m 2 ~ e 4 T for n — 1, 
where M = 4 supersymmetry is restored in the e — > limit. 

Although perturbation theory seems to have failed, we find it plausible that the heavy 
fields should play a role in setting the scale of supersymmetry breaking at y/eT. This 
scale is just the geometric mean of the tree- level masses of the light and heavy fields, 
but its relevance to the dynamics may emerge in field theory only through features 
peculiar to strong coupling gauge interactions. 

To summarize our discussion: In no sense should spinning M5-brane solutions with 
a large angular momentum be regarded as the analog of improved lattice actions for 
pure QCD 4 . The fermions' tree level masses are not all far above the confinement scale. 
In fact, some of them are far smaller than this scale. We do not understand the field 
theory mechanism that sets the scale of supersymmetry breaking and confinement at 
0(y/eT) for the case of one large angular momentum. However, if we take this scaling 
as a given, we do understand why glueball mass ratios should be universal for two large 
unequal angular momenta and somewhat different when they are equal: it is because 
the pattern of supersymmetry breaking is universal in the first case, but different in 
the second. It would be delightful if the near-restoration of supersymmetry could 
be exploited to make more definite predictions about the confining dynamics than is 
possible for pure QCD 4 . 

The D3-brane with periodic Euclidean time has been proposed as a model for pure 
QCD 3 [23]. Without angular momentum, the bare three-dimensional fermion masses 
and also the mass gap are on the order of T. With n equal large Euclidean angular 
momentum, the light fermion masses are mbarc ~ cT, whereas the mass gap is M gap ~ 
e 1_ iT for n — 1 or 2. Again it is not easy to see what in field theory sets this 
confinement scale. Universality of glueball mass ratios can be understood in a similar 
way to the four-dimensional case, but we regard it as a more pressing problem to 
understand the overall scaling of the mass gap from a field theory perspective. 

The supergravity analysis of the previous subsection applies to Euclidean spinning 
M2-branes as well. However, in view of the poor understanding of the large N field 
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theory we refrain from drawing any further comparisons. 

V. CONCLUSION 

The purpose of this paper was to provide a comprehensive study of the thermo- 
dynamics of spinning D3-, M5-, and M2-branes. Our results address the following 
aspects: 

• Thermodynamic Stability and Critical Behavior 

We derived general stability constraints for near-extreme D3-, M5- and M2- 
branes, all with the maximum number of independent spins. The region where the 
entropy is subadditive as a function of both the energy and the angular momenta 
is strictly smaller than the region where the specific heat is positive. To illustrate 
one of the possible instabilities that develops when the entropy is not subaddi- 
tive, we studied states where the angular momentum is unevenly distributed over 
the world volume, and found them to be entropically preferred over the uniform 
distribution despite the fact that they are not themselves stable configurations. 

• Field Theory 

We generalized the field theory model of [1] to the case of multiple R-charges; 
the dynamics is proposed to be described effectively by a collection of weakly 
interacting massless "quasi-particles" whose R-charges and spins are those of the 
abelian theory on a single brane. If one assumes in a addition that there is a 
nonzero "permeability" for applied "voltages" dual to the various R-charges, then 
this model reproduces the critical behavior of the supergravity analysis. 

• Euclidean spinning branes 

Euclidean branes with at least one angular momentum density which is large 
compared to the energy density restore some fraction of supersymmetry. We 
showed how this supersymmetry restoration arises in the field theory: an element 
of the R-symmetry group is inserted in the partition function, which for one large 
spin approaches (— 1) F . The spectrum is nearly super symmetric, with the Bose- 
Fermi degeneracy broken by AM ~ eT, where e is a dimensionless parameter 
that goes to zero in the limit of large angular momentum. The supergravity 
analysis reveals that if there is confinement, it occurs on a scale M gap ~ j^ 1 -™/ 2 
when n equal Euclidean angular momenta are made large. We have explained 
how this picture elucidates some of the results of previous papers [6,11]. We have 
left open the problem of identifying the field theory mechanism that sets the scale 
of supersymmetry breaking. 
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